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^ ■ ABSTRACT 

bJOl Context. Multi-wavelength time-series observations with high cadence and long duration are needed to resolve and understand the 

^ . many variations of magnetically active late-type stars, which is an approach often used to observe the Sun. 

<^ ■ Aims. We present a first and detailed study of the bright and active KOIV-III star HD 123351. 

Methods. We acquired a total of 955 high-resolution STELLA echelle spectra during the years 2006-2010 and a total of 2,260 

CO \ photometric Vic data points during 1998-2010. These data are complemented by some spectra from CFHT and KPNO. 

■ Results. The star is found to be a single-lined spectroscopic binary with a period of 147.8919+0.0003 days and a large eccentricity 
of ^=0.8086+0.0001. The rms of the orbital solution is just 47 ms~\ making it the most precise orbit ever obtained for an active 
binary system. The rotation period is constrained from long-term photometry to be 58.32+0.01 days. It shows that HD 123351 is 
a very asynchronous rotator, rotating five times slower than the expected pseudo- synchronous value. Two spotted regions persisted 

. throughout the 12 years of our observations. We interpret them as active longitudes on a differentially rotating surface with a AP/P 

of 0.076. Four years of Ha, Can H&K and Hei D3 monitoring identifies the same main periodicity as the photometry but dynamic 
spectra also indicate that there is an intermittent dependence on the orbital period, in particular for Can H&K in 2008. Line-profile 
inversions of a pair of Zeeman sensitive/insensitive iron lines yield an average surface magnetic-flux density of 542+72 G. The time 
Q , series for 2008 is modulated by the stellar rotation as well as the orbital motion, such that the magnetic flux is generally weaker 

■ during times of periastron and that the chromospheric emissions vary in anti-phase with the magnetic flux. We also identify a broad 
^ ' and asymmetric lithium line profile and measure an abundance of log;i(Li) = 1.70 + 0.05. The star's position in the H-R diagram 
^ \ indicates a mass of 1.2+0.1 M© and an age of 6-7 Gyr. 

Conclusions. We interpret the anti-phase relation of the magnetic flux with the chromospheric emissions as evidence that there are 
two magnetic fields present at the same time, a localized surface magnetic field associated with spots and a global field that is oriented 
towards the (low-mass) secondary component. We suggest that the inter-binary field is responsible for the magnetic-flux dilution at 
^ ' periastron. It is also likely to be responsible for the unexpected slow and asynchronous rotation of the primary star. 

: 

_ Key words, stars: activity of - starspots - stars: individual: HD123351 - stars: magnetic field - stars: late-type 

QQ 1. Introduction high eccentricity, may allow us to prove or disprove the specula- 

^ t . .... , . „ tion that a dynamo would operate more efficiently at periastron 

^ Cool main-sequence stars in binaries appear to be magnetically ^^^^ ^^^^^ ^j^^^^^^ ^^^1^ ^^^^ ^ ^ 

^ very active when the orbital period is as short as a few days, influenced by tidal forces (e.g. Schrijver & Zwaan 1991). Such 

' Evolved components in binaries are suspected to display photo- ^ ^. ^ presented and analyzed in the present paper. 

, 1^ spheric spots and strong chromospheric emission at significantly rj^^ ^^^^ J2335 1 (a = 14''06'"26^ 6 = 30°50'47" [2000] 

a longer periods, say, 20 days. This has been generally interpreted y ^ was found to exhibit strong Ca 11 H&K emission as 
as meaning that there is a stellar dynamo that needs a convective Vienna-KPNO Doppler imaging candidate survey 

_Cd envelope with a long^enough turn-over time and short-enough (strassmeier et al. 2000). The survey spectra allowed us to deter- 
rotation to function efficiently (Durney & Latour J978j. The ex- emission-line fluxes, a projected rota- 

istence of active longitudes in tidally interacting binaries is part ^j^^^^j ^^j^^j 5 2±2 km s'^ , a logarithmic lithium abundance 
of the proof of this concept (see, e.g., Berdyugina 2007). The ^ 9^ .^^^ ^^^j^j ^^^^^.^^^ ^ ^^^^ 

observational evidence is less clear for stars in a binary where ^-^^^^ ^ ^^^ ^-^^^^ ^^^^^ ^^^^^ ^j^^^ ^^^^ ■ ^-^^ ^^^^ 

the rotation and the orbital revolution is not synchronized or, , ^ , ^u.^^ « .^^^ ^ i^^^^ u^^^ 

, . „ , . I t . had observed a single star. To our knowledge there had been no 

u^^^^'n'^nnTPc having additionally a highly eccentric or- ^^.^^^^^ ^j^^ ^^^^^^^^^ ^^^^ ^j^^ ^ spectroscopic or a vi- 

bit (Hall 1991 , Strassmeier 2001 ). A binary system with an ac- ^^^j ^. ^ -^^ -^^ ^^^^^^^ brightness. Our new data in 

tive component and all of above properties, in particular wiffi a ^j^^ p^^^^^^^ p^p^^ -^^^^^^ ^ single-lined spectroscopic 

Send offprint requests to: K. G. Strassmeier binary with a very eccentric orbit. 

* Based on data obtained with the STELLA robotic telescope in The light variability appears to have been first detected by 

Tenerife, an AIP facility jointly operated by AIP and lAC, and the the Hipparcos satellite which led Kazarovets et al. (1999) to in- 

Potsdam Automatic Photoelectric Telescopes (APT) in Arizona, jointly elude HD 123351 (HIP 68904) in the 74th name list of variable 

operated by AIP and Fairborn Observatory. stars as CZ CVn. Simbad identifies it as a "pulsating star" with 
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evidence neither for nor against it. The Hipparcos Hght curve 
was analyzed a few years later by Koen & Eyer f2002), who 
found a photometric period of 59.4 days with an amplitude of 
CP028. A year earlier we had published our first three years of 
Automatic Photoelectric Telescope (APT) VI data from 1998- 
2000 (in a paper on APT data reduction by Granzer et al. 12001b 
that also clearly showed the 60-day periodicity. 

The star was detected as an X-ray and EUV source by EUVE 
and ROSAT with count rates of 0.88 and 0.021, respectively 
(Lampton et al. J997 ). It is Hsted in the ROSAT all-sky sur- 
vey (Schwope et al.^000) and identified in the Hamburg/RASS 
bright source identification catalogue (Zickgraf et al. 120031) . A 
study of the X-ray selected Seyfert galaxies by Vaughan et al. 
(120011) Hsted HD 12335 1 in their Galactic sample with a Galactic 
column density of 1.25x10^^ cm"^. A study of the local kine- 
matics of K and M giants with Coravel/Hipparcos/Tycho data 
(Famaey et al. 2005 ) included HD 123351 but did not provide 
a new radial-velocity measurement. The new Pulkova radial- 
velocity catalog of HIP stars lists HD 123351 with a single mea- 
surement of -1 1.2 ± 0.8 km s"^ (Gontcharov 2006). 

Therefore, the star appeared to be an interesting and rather 
bright chromospherically active star and deserved a more de- 
tailed study. In this paper, we present and analyze an extensive 
time- series database consisting of high-resolution echelle spec- 
tra obtained as part of the science demonstration of the first of 
our two robotic STELLA telescopes in Tenerife, complemented 
by some previous KPNO and CFHT spectra, and 12 years of 
photometric monitoring with our APT in southern Arizona. In 
Sect. 121 we describe our data. In Sect. 3, we focus on the recovery 
of the star's basic astrophysical parameters such as the rotational 
period, the orbital elements, stellar mass, radius and age, lithium 
abundance, Rossby number, and evolutionary status in general. 
Section [4] presents the time-series analysis of the photometric 
and spectroscopic data in order to reconstruct the stellar surface 
structure as a function of rotational and orbital phase. Section [5] 
is an attempt to gather the various results together. Our summary 
and conclusions are presented in Sect. [6l 



Table 1. Spectroscopic and photometric observing logs. 



Observatory 


Type 


Instr. 


Time range 










JD 245+ 




STELLA-I 


Echelle spec. 


SES 


3912-5213 


955 


KPNO CF 


Coude spec. 


Coude 


0908-1673 


20 


CFHT 


Coude spec. 


Gecko 


1683-1687 


4 


APT T7 


Vic 


PMT 


0911-5200 


2260 



Note. N is the total number of individual spectra or light-curve points. 
CF denotes the 0.9m Coude Feed telescope at Kitt Peak. PMT means 
PhotoMultiplier Tube. 



of 2260 V and Ic data points were obtained until the end of 
the observing season 2010 and are presented in this paper in 
Fig. [U The integration time was 10 sec in V and Iq- All mea- 
surements were done diff'erentially with respect to HD 123929 
(G8, y=7^ 27, B - y-0.85, Oja 1985; ¥^=7^360, B - y=0^835, 
H0g et al. I2000b . The check star was HD 123533 (KO). Both 
stars appear constant to within the telescope's internal precision. 
The APT data were transformed to absolute values with nightly 
observations of up to 60 standard stars. For further details of the 
observing procedure and the APT data reduction in general, we 
refer to Granzer et al. (120011) . 

2.2. STELLA/SES spectroscopy in 2006-2010 

High-resolution time- series spectroscopy was obtained with 
the STELLA Echelle Spectrograph (SES) at the robotic 1.2-m 
STELLA-I telescope in Tenerife, Spain (Strassmeier et al. 120041 
[20 10a). A total of 955 echelle spectra were acquired over the 
course of four years. The SES is a fiber-fed white-pupil echelle 
spectrograph with a fixed wavelength format of 388-882 nm. 
Despite increasing inter-order gaps in the red, it records the 
range 390-720nm continuously. The instrument is located in a 
dedicated room on a stabilized optical bench and is fed by a 12- 
m long 50yum Ceram-Optec fibre, corresponding to an entrance 
aperture of 1.7'' on the sky. This fiber enables a two-pixel res- 
olution of 7?=55,000. The spectrograph's heart is a 31 lines per 
mm R2 grating. Two off'-axis parabolic collimators, one fold- 
ing mirror, and a prism as a cross disperser transport the light 
into the f/2.7 kata-dioptric camera with a 20cm corrector and 
a 40cm spherical mirror. The CCD is currently an E2V 42-40 
2048x2048 13.5yum-pixel device. Its quantum eflftciency is 90% 
at 650nm, and 65% at 400 and 800nm, respectively. Together 
with a second-generation CUO controller, the read-out noise is 
3-4 electrons rms at 200 kbyte/s. A closed cycle cooler keeps 
the detector cooled to -130°C. 

Integrations on HD 123351 were set to exposure times of 
between 3600 sec and 1200 sec and achieved a signal-to-noise 
ratio (S/N) of between 200:1 to 50:1 per resolution element, re- 
spectively, depending on weather conditions. One spectrum on 
almost every clear night was obtained between June 25, 2006 
and June 22, 2010. Numerous radial velocity standards and stel- 
lar comparison targets were also observed with the same set- 
up. Earlier data in 2006 and during March- October 2007 had to 
be shifted by +280 m s"^ because of a fiber misalignment prob- 
lem. The time-series velocities are plotted in Fig. O The best 
rms radial- velocity precision over the three years of observation 
was 30 m s"^ . The time series also included nightly flat-field ex- 
posures, bias readouts, and Th-Ar comparison-lamp exposures. 
For details of the echelle data reduction, we refer to Ritter & 
Washuettl (12004b and Weber et al. ( 2008 ) and our previous paper 
on HD 1 (Strassmeier et al. 1201 Oct . The data log is summarized 
in Tablelll 



2. Observations 

2.1. APT photometry since 1998 

High-precision photometry of HD 123351 was obtained with 
the former University of Vienna 0.75m twin automatic photo- 
electric telescopes (APTs) at Fairborn Observatory in southern 
Arizona (Strassmeier et al. 1 1997 1 Granzer et al. 2001). The T7- 
APT Amadeus was used since JD 2,450,900. It achieved an ex- 
ternal precision of 3.2mmag in Johnson-Cousins Vic- A total 



2.3. KPNO/CF spectroscopy in 1998 and 2000 

Twenty KPNO observations were obtained with the 0.9-m coude 
feed (CF) telescope in April 1998 and in March 2000 (Table[B. 
The F3KB 3k-CCD was used together with grating A, the long 
collimator and a 280 jim slit. This configuration allowed a re- 
solving power of R ^^28,000 at a dispersion of 4.81 km s"Vpixel. 
The wavelength coverage was 30 nm and the spectra were cen- 
tered at 649 nm. The exposure time was 1200 seconds resulting 
in an average S/N of 100:1. 
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1998 



1999 



2000 



2001 



2002 



2003 



2004 



2005 



2006 



2007 



2008 



2009 



2010 



7.30 



7.40 



7.50 



7.60 



7.70 



7.80 



1.00 
1.02 
1.04 
1.06 
1.08 
1.10 
1.12 
1.14 
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Fig. 1. The long-term brightness variations of HD 123351. The top panel shows our twelve years of V-band photometry (dots), the 
lower panel plots the corresponding V - I color index. The lines are the seasonal fits of the two-spot models from Sect. IH These 
models are plotted here as spherical projections at four diff'erent phases (from top to bottom, 0°, 90°, 180°, 270°). The gray scale 
indicates the spot temperature. Note the rotational modulation with a period of ^^60 d and a long-term variation in both V and V - 1 
in the sense that the star appears redder when fainter. 



All KPNO spectra were reduced and extracted using the 
Interactive Reduction and Analysis Facility (IRAF) provided 
by NOAO. More details of the standard coude- spectra reduc- 
tion procedure were given in previous papers, e.g. by Weber & 
Strassmeier (119981) . 

2.4. CFHT/Gecko spectroscopy in 2000 

Three very high-resolution lithium spectra were obtained con- 
secutively on one night at CFHT with Gecko and an order- 
sorting filter in the 9th order at HJD 2,451,687.75. A fourth inte- 
gration of 900 sec was taken centered at 644 nm at 2,45 1 ,683.73. 
Equipped with the CAFE fiber feed module Gecko provided a 
spectral resolution of 120,000. CAFE feeds light through a 27m- 
long lOOyum Ceram-Optec fiber to a four-slice Bowen-Walraven 
image sheer before it enters the collimator. In combination with 
the 4600x2048 13.5yum-pixel EEVl CCD, it gives a lOnm wave- 
length coverage centered at 671 nm. An integration time of 
600 sec per exposure gave a peak S/N of 300:1 per resolution 
element. For a brief description of the data reduction, we refer to 
the paper by Strassmeier & Rice (I2003i) . The data log summary 
is given in Table [T] 



3. Stellar parameters of HD 123351 

3.1. Orbital elements 

Radial velocities of the STELLA spectra were determined from 
an order-by-order cross correlation with a synthetic G8 template 
computed from a Kurucz (il993i) atmosphere. Only the central 



19 echelle orders are used, providing a wavelength range of 
^^150nm. The coude spectra from KPNO and CFHT were an- 
alyzed similarly but using the IRAF cross-correlation program 
fxcor together with various reference- star spectra. All veloci- 
ties are corrected for the solar- system's barycenter (the full data 
are given in machine readable format in the electronic attach- 
ment). The gravitational redshift is expected to be 133 m s"^ ac- 
cording to the mass and radius in Table [3] but was not corrected 
for. 

For the orbital solution with non-zero eccentricity we use 
the prescription from Danby & Burkardt (119831) to calculate 
the eccentric anomaly. The elements are solved for by using 
the general least- square fitting algorithm MP FIT (Markwardt 
2009 ). More details about the orbit computation are discussed in 
Weber & Strassmeier (2011) and we refer the reader to this pa- 
per. Owing to unclear zero-point diff'erences between the three 
data sets, and the diff'erent set ups, we include the earlier ve- 
locities from KPNO and CFHT only for the period determina- 
tion. An SBl solution of all velocities, but with zero weight 
given to the KPNO and CFHT data, resulted in an eccentricity of 
0.809+0.001 and a period of 147.890 d. We note that all previ- 
ous KPNO and CFHT velocities fell into the phase range 0.1-0.8 
where almost no velocity variations because of the high eccen- 
tricity. From these spectra alone, one could conclude that the star 
has a constant radial velocity, as actually done by Strassmeier et 
al. (|2000). For the initial orbital solution, a total of 955 STELLA 
radial velocities were used. The standard error in an observation 
of unit weight was 73 m s"\ which is the best so far achieved for 
an orbit of a chromospherically active star (e.g. Fekel et al. 119991 
Griffin |2009» Torres et al. i2009} . However, this rms includes the 
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Phase JD 

Fig. 2. Radial velocities of HD 123351. a. Phased STELLA radial velocities with the orbital period of 147.8919 d. Dots are the 
observations and the line is the orbital fit to the spot-corrected elements in Table [21 b. Top panel: The four years of monitoring. 
Such data sampling is made possible with robotic telescopes. The line is again the orbital solution from Table [2] (column "Spot 
corrected"). Lower panel: Residuals after removal of the orbital fit. The modulation is due to stellar rotation and the insert shows a 
DFT periodogram with the strongest peak at 61.0+1.1 d, in agreement with the photometric period. The full amplitude varied over 
the years between 400 m s"^ in 2007 and 100 m s"^ in 2009, indicating the variable asymmetry of the spot distribution. The line is 
a Blackman-function fit to the data. 



modulation due to the asynchronously rotating spotted surface. 
Table [2l (left column) lists the orbital elements from the unal- 
tered initial data. In Sect. 13.21 we refine the orbit by removing 
the star spot jitter from the radial velocity data and then achieve 
a standard error of 47 m s"^ close to the measuring error of ap- 
proximately 30 m s"^ Fig. [2^ plots the spot-corrected radial ve- 
locities along with the orbital fit, Fig. [2]3 (top panel) shows the 
unaltered original data. The diff'erence cannot be seen by eye. 

We note that the low mass function, together with our in- 
ability to see a secondary star in our spectra, suggests that the 
orbit has a low inclinationwith respect to the sky. If we assume 
at least a 2^5 brightness diff'erence between the secondary and 
the primary at red and blue wavelengths, the secondary can- 
not be of earlier spectral type than G2(V). Adopting its mass as 
m2=l.O Mq, the observed /(m)=0.01 13 provides a lower limit to 
the inclination of almost 22°. A more likely secondary would be 
a red dwarf of mass m2 ~ 0.5 M©, which suggests a more likely 
inclination of ^ 40° and a minimum separation of ^ 2 10^ km 
(^^5 stellar radii) during periastron. 

3.2. Rotational parameters 

The photometric data have high amplitude, both short-term and 
long-term, variations. The short-term variations are clearly due 
to rotational modulation, while the long-term variation is likely 
due to some sort of multi-periodic spot cycle. 

A photometric period of ^^60 d is obtained with a very high 
level of confidence, which we interpret to be close to the true 
stellar rotation period. This value is based on a one-dimensional 



Table 2. Spectroscopic orbital elements. 



Parameter 


Original 


Spot corrected 


Porb (days) 


147.890+0.001 


147.8919+0.0003 


T'periastron (HJD245+) 


3980.616+0.004 


3980.6098+0.0021 


y (km s~^) 


-8.703+0.003 


-8.7057+0.0017 


Ki (km s-^ 


15.35 +0.03 


15.3518+0.0091 


e 


0.809+0.001 


0.8086+0.0001 


CO 


15.34 +0.01 


15.326+0.018 


ai sin / (10^ km) 


18.35 +0.03 


18.368+0.012 


/(m) (Mo) 


0.01128+0.00010 


0.011317+0.000023 


fit rms (m s ^ ) 


73 


47 



Note. The columns are for the unaltered data ("Original") and for 
star spot-jitter removed data ("Spot corrected"). 



discrete Fourier transform (DFT) analysis of the V/c-band data. 
The best-fit twelve-year average period is 61.9+0.4 d from a 
Gaussian fit to the V and /c-band DFT peaks (the results from 
y data are shown in Fig. [3^). However, other significant peaks 
show up at slightly higher frequencies even in the pre- whitened 
data. Repeating the analysis for three data subsets by splitting the 
entire time series into three parts of similar length, we get peri- 
ods of 61.4 d for the first (1998-2001), 62.0 and 57.3 d for the 
second (2002-2005), and 59.2 d for the third (2006-2010) part 
of the data set, identifying the variable and sometimes even two 
significant periods. The typical error is 0.7 d. For spotted stars, 
this is a sign of surface diff'erential rotation, i.e., spots modu- 
late the light with diff'erent periods when they appear at diff'erent 
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latitudes on the stellar surface. It also means that we cannot sim- 
ply interpret the photometric period to be the rotational period 
of the star. Seasonal changes occur when spots disappear after 
a certain (life)time but others appear, possibly cyclically as on 
the Sun, and possibly at different latitudes. If either the cycle is 
short enough or the database of sufficiently long duration and 
the sampling dense enough, it allows us in principle to deter- 
mine the amount of differential rotation (e.g. Olah et al. l2QQ3b . 
Its sign, however, could only be determined from simultaneous 
Doppler imaging. 

To confirm the above-mentioned period splitting, and pin- 
point its onset and duration, we performed a two-dimensional 
time-frequency analysis with the program package Time- 
Frequency Analysis (TiFrAn) (Kollath & Csubry 2006). The 
time-frequency analysis is conventionally applied to continuous 
datasets, whereas our observations are gapped. Thus, the yearly 
interruptions are filled with straight lines, which do not alter our 
main results but may produce spurious features. More details 
about the procedure can be found in Kollath & Olah ( 2009b . We 
chose the short-term Fourier transform (STFT) but also com- 
pared the results with those obtained using the Choi- Williams 
and the pseudo-Wigner kernel. The resulting patterns of peri- 
ods are identical to those of the above DFT analysis to a very 
high degree. Fig.fSj) shows the frequency distribution over the 12 
years of data from the STFT analysis. The continuous variability, 
and the splitting of the rotational period, is clearly visible. The 
splitting started in 2003 and lasted until 2005. The 2004 observ- 
ing season is the one where the light curve shape became double- 
humped and the V and Ic amplitudes reached their smallest value 
in our entire time series. At the same time, the photometry ap- 
peared bluest, indicating, together with the low amplitudes, that 
the spot coverage had reached a minimum. 

The coadded R = 120,000 CFHT spectrum with S/N=500 
in the wavelength range 63 8-646 nm was measured for spectral- 
line broadening. We followed the recipe of Fekel (1997) and 
found a best-fit vsin/ value of 1.8+0.7 kms"^ adopting a typ- 
ical radial-tangential macroturbulence for a K subgiant/giant of 
3 km s"^ We note that such a low v sin / cannot be resolved by 
our spectra and thus formally constitutes only an upper limit. 
Together with Pphtm ~ ^rot ~ 60 d and the assumption that the 
rotational axis is perpendicular to the orbital plane, we infer that 
the minimum stellar radius is R sin / ^^2. 1+0.8 R©. 

3.3. Rotation induced radial-velocity variations 

Today's high-precision radial- velocity surveys of stars with ex- 
oplanets may record the rotation of the host star if it has dark 
spots just large enough to mimic a "transit" (Saar & Donahue 
11997) . Such rotational modulation was first detected in active 
pre-main sequence stars with large cool spots (e.g. Stelzer et 
al. 12003 ^) but also in an active main- sequence star (Ruber et al. 
120091) . Amplitudes of 2 kms"^ were detected for the extremely 
active T Tauri star LkCa 19 (Huerta et al. 120081) , while a flare on 
the M-dwarf CN Leo introduced radial velocity jitter of several 
tens of ms"^ (Reiners r2009b . We found radial velocity varia- 
tions with an amplitude of 270 ms"^ for the single, active giant 
3 1 Com with a period in agreement with the photometric period 
from MOST satellite data (Strassmeier et al. i2010bb . 

HD 123351 exhibits radial- velocity modulations that are 
even visible by eye in the plot of the residuals in Fig.[2]3 after the 
removal of the orbital variation. The peak-to-peak amplitude var- 
ied across the range 300-400 m s"^ in 20 07 and was 100 m s"^ 
in early 2009. The simulations of Hatzes (120021) suggested that 
spots induce a radial- velocity amplitude of at least hundred m/s 
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b) STFT time-frequency analysis 

Year 

1998.51 2001.25 2003.98 



2006.72 



0.004 



2009.46 



f [c/d] 
0.02 

0.015 

0.01 

0.005 




2451000.0 2452000.0 



2453000.0 
JD 



days 
50.0 

66.6 

100. 

200. 

inf. 



2454000.0 2455000.0 



Fig. 3. Periodograms from HD 123351 photometry, a. Discrete 
Fourier transform (DFT) periodogram from the pre- whitened V- 
band data (top panel). The twelve-year average period is 61.9 d 
(0.01616 c/d). The bottom panel is the corresponding window 
function. The insert in the top panel shows the entire search 
range, b. Two-dimensional short-term Fourier transform (STFT) 
analysis. The power distribution in the time-frequency plane is 
shown as a gray scale (black means highest power, white no 
power). It shows that the photometric period reached values be- 
tween a maximum of 62.0 d and a minimum of 57.3 d. Around 
2,453,000 (2002-2005), a bimodal distribution was recovered. 
Note that the excess power in the lower left corner is due to the 
length of the data vector and is artificial. 



for even low-activity stars. His relationship of velocity amplitude 
Sv in m s"^ with spot coverage A in per cent and line broadening 
V sin / in km s"^ (^v = (8.6 v sin /- 1.6)A^-^) would suggest a vari- 
able spot coverage in the range of A 9 - 42 % for HD 123351. 
This is incorrect by a factor of three when compared to the more 
detailed simulations by Desort et al. (2007), who found 3- 12 % , 
mostly because these authors took into account the spot tempera- 
ture. However, this diff'erence is of no relevance here because we 
only wish to make an order-of-magnitude estimate. Makarov et 
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Table 3. Summary of astrophysical data of HD 123351. 



Parameter 



Value 



Maximum V magnitude 
Spectral type 
Hipparcos distance 

My 

Av 
V sin / 

Rotation period 

Inclination 

Radius 

Luminosity 

Mass 

Age 

[Fe/H] 

Lithium abundance 
Average turn-over time 
Microturbulence 
Macroturbulence 
Average magnetic flux 



7^235 
KO IV-III 
97.8^^ PC 
+2?20+0.14 
(r08 

4780+70 K 
3.25+0.30 
1.8+0.7 km s-i 
58.32+0.01 d 

^40° 

15.4:2-8 Lo 
1.2+0.1 Mo 
^6-7 Gyr 
0.00 + 0.08 
1.70+0.05 
40 d 

1.15kms-i 
3.9 km s"^ 
542+72 G 



al. (2009) extended these analytical relations to astrometry. We 
performed a DFT analysis of the radial- velocity residuals shown 
in Fig. [2I3 and found the best-fit period to be 61.0+1.1 days, in 
agreement with the photometric period. 

The rotation of the spotted stellar surface introduces system- 
atic jitter into the radial velocity curve. However, this jitter is 
unevenly distributed across one orbital revolution because of the 
asynchronous stellar rotation. We note that the residual radial- 
velocity curve appeared even quintuple humped in 2008, just as 
the light curve appeared double humped at the same time. We 
remove the starspot jitter from the original data and then repeat 
the determination of the orbital elements from Sect. 13.11 This is 
achieved by convolving the residuals in Fig. [2J3 with a smooth- 
ing function (Blackman function) and then subtracting it from 
the radial- velocity data. The largest values that were added to 
or removed from the data were -hlOO m s"^ and -200 m s"\ re- 
spectively. The revised orbital elements are given in the column 
"Spot corrected" in Tabled 

The eff'ect of the applied procedure on the orbital elements 
is also demonstrated in Table O We note that the elements did 
not change much, i.e. the corrected solution is largely within the 
errors of the uncorrected solution. Typically, the error per ele- 
ment became smaller by a factor 2-3. This is what is expected if 
the jitter due to rotational modulation spreads equally, or nearly 
equally, around a mean radial velocity. For our long-term obser- 
vations, this is indeed the case because the phased orbital radial- 
velocity curve is sampled in total by 25 stellar rotations (for 
each orbital revolution we have 2.5 stellar rotations). Because 
the starspot distribution changes with time, we do not just cor- 
rect with a harmonic sinusoidal of period Prot and fixed ampli- 
tude but with the actual residuals of the rotationally modulated 
radial velocity. 

3.4. Luminosity, radius, metallicity, mass, and age 

Fig. |4] compares the position of HD 123351 in the H-R diagram 
with theoretical evolutionary tracks. We note that the revised 
Hipparcos parallax of 10.22+0.67 mas (van Leeuwen i2007j cor- 
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Fig. 4. The position of HD 123351 (dot) in the theoretical H- 
R diagram (logL/L© vs. logT in K). Shown are evolutionary 
tracks for masses of 1.1-2.5 M©, solar metallicity, and core 
overshooting (dotted lines; Pietrinferni et al. 2004). The best- 
fit model track (full line) suggests a mass for HD 123351 of 
1.2+0.1 Mq. The insert is a zoomed image of the box around 
the observed location indicating the nominal ages of the 1.2-Mo 
track (again the full line) and a l.l-M©, 1.3-Mo, and 1.5-Mo 
track (dotted lines), respectively. This suggests a most-likely age 
ofHD 123351 of 6-7 Gyr. 



responds now to a distance of 97. 8 ^ pc, which is diff'erent by 
-hl.7 pc from the original data reduction. At this distance, it is 
only moderately aff'ected by interstellar extinction. We adopted 
the mean extinction value from Henry et al. (12000b for the abso- 
lute magnitude calculation (Ay =0.8 mag kpc"^ and E(B - V) = 
Av/3.3). The brightest recorded V magnitude of 7^235 con- 
verts into an absolute visual magnitude of My = -h2?20 + 0.14. 
Its dereddened B -V based on the Tycho color (1^029+0.01) 
is 1^00. This can be converted to a T^ff of ^4840 K with the 
transformation of Flower (1996). The bolometric magnitude of 
HD 123351 is +1^78 obtained with a solar bolometric mag- 
nitude of+4^75 and a bolometric correction (B.C.) of -0.42 
(Flower [1996 ). The +70 K uncertainty in Tqq (see next para- 
graph) corresponds to an uncertainty in the B.C. of +0.04 and 
the star's luminosity becomes 15.4^2 3 L©. These errors are still 
driven by the error in the parallax rather than in the eff'ective 
temperature (+70 K). The Stefan-Boltzmann law gives a stel- 
lar radius of 5. 7 4'^^-^^ R©, which is appropriate for a KO sub- 
giant/giant. The minimum radius inferred from v sin / and Prot in 
Sect. 13.21 of 2.1+0.8 R© suggests a minimum inclination of as 
low as / ^ 23°. 

Our SES spectra are used to refine the eff'ective tempera- 
ture, gravity, and metallicity. We employ the spectrum- synthesis 
technique and the numerical tools described in Allende-Prieto et 
al. ( 2004) with a grid of synthetic ATLAS -9 spectra tailored to 
the stellar parameters of HD 123351. Five spectral orders around 
600 nm were chosen and the result per spectral order combined 
on the basis of a weighted least- squares minimization. The av- 
erage and the standard deviations then constitute our final val- 
ues and their internal precisions. We found Tq^ = 4780 + 70 K, 
logg=3.25+0.30, and [Fe/H]=0.00 + 0.08 with a prefixed value 
for the microturbulence of 2 kms"^ and a total line broadening 
of 3.80 km s-^ 

These values are compared with theoretical evolutionary 
tracks in Fig.|4l We selected post main-sequence tracks from the 
grid of Pietrinferni et al. (12004b . All tracks are for solar metallic- 



6 



K. G. Strassmeier et al.: Binary-induced magnetic activity? 



ity and core overshooting. The HD 123351 position suggests a 
most Hkely mass of 1.2+0.1 M©. We note that the models with- 
out core overshooting would increase the mass by about 0.1 M© 
and decrease the age by ^ 1 Gyr. The star appears on the ascent of 
the red-giant branch (RGB) at an age of approximately 6-7 Gyr. 
The ^^C^'^N lines around 800.35 nm are clearly detectable with 
a combined equivalent width of 96 mA. However, the equiva- 
lent width of the ^^C^'^N line at 800.46 nm is at ^2.5 mA very 
weak and probably below our detection limit but provides a 
lower limit to ^^C/^^C of ^40, in agreement with the expectation 
that small ^^C/^^C ratios (around 10-15) are only found in giants 
that have evolved up to or beyond the RGB (e.g. Allende Prieto 
et al. [2004"). At this point we recall the findings of radius and 
temperature inconsistencies for active low-mass eclipsing bina- 
ries with model computations (e.g. Morales et al. |2009l 12010 ). 
Oversized stars with a temperature deficit caused by chromo- 
spheric activity could alter the true stellar position in the H-R 
diagram, thus aff'ect the mass and age estimate. 

3.5. Lithium abundance 

We employ the averaged 7?= 120,000 CFHT spectrum for the 
measurement of the Li equivalent width (Fig. [5]). This spec- 
trum has a peak S/N of almost 500:1. A Gaussian fit gives an 
equivalent width of the combined 670.78 nm line (^Li+^Li) of 
90+3 mA (the error is an estimate inferred from the spread in 
the values from the three individual spectra and the quality of 
the Gaussian fit). Subtraction of a reference spectrum with pre- 
sumably no lithium (16 Vir in Fig. O removes a contribution 
of 16 mA from the nearby Fe i+V i blend. The remaining 74 mA 
converts into an abundance of log ^(Li)=1.70 (log^(H)=12) with 
Teff = 4780 K and log g = 3.0 from the NLTE tables in Pavlenko 
& Magazzu (119961) . Its uncertainties produced by the measuring 
error and the error of +70 K in Tq^ and +0.3 in logg, combine 
to ^0.05 dex. 

The very low rotational broadening of v sin/ = 1.8 kms"^ 
enables a detailed line-profile modeling of the ^Li/^Li iso- 
tope ratio. As this is a complex synthesis based on a three- 
dimensional, radiative hydrodynamics simulation, we refer to 
a forthcoming paper by Steff'en et al. (120111) . Not much is 
known about this ratio in magnetically active stars (see, e.g., 
Balachandran et al. 12000b but it plays an important role in con- 
straining or even excluding some mixing mechanisms and deter- 
mining the likelihood of extra lithium production. HD 123351 
shows a moderately strong but broader than expected Li i line 
at 670.8 nm (see Fig.O. The line strength is weaker than to the 
nearby Fe i lines at 67 1 .03 nm and 670.5 1 nm, i.e. 0.76 compared 
to 0.68 and 0.67, respectively, but the equivalent width is larger 
by more than 20%. A weak excited Li feature at 610.36 nm is 
also present but strongly blended with the nearby double Fei 
line and appears rather uncertain for a quantitative analysis. We 
measure an equivalent width of 8+4 mA from a double Gaussian 
fit, close to the detection limit of, say, 3-5 mA. Therefore, we re- 
frain from modelling this line. 

3.6. Rossby number and dynamo activity 

The Rossby number, Ro = Pm/^c, is widely used as an indicator 
of the existence and efficiency of a dynamo process. While the 
rotation period can be observed very accurately, the convective 
turn-overtime, r^, remains an elusive parameter. It describes the 
ratio of a characteristic length of convection to a characteristic 
velocity, ^/w, and is either determined theoretically (e.g. Kim & 



HD123351 (HJD2,451,687.75) 




6702 6704 6706 6708 6710 6712 

Fig. 5. Lithium 670.8-nm region in the spectrum of HD 123351. 
Shown are three consecutive R = 120,000 spectra obtained at 
CFHT in May 2000 at orbital phase 0!496 along with a refer- 
ence spectrum of the inactive K0.5III star 16 Vir. All spectra are 
shifted in intensity for a clearer display. A very high-resolution 
spectrum of the Sun from the National Solar Observatory solar 
atlas (Kurucz et al. 11984b is shown for comparison. HD 123351 
shows a significant lithium line with a total equivalent width of 
90 mA. 



Demarque l 1996b or semi-empirically (e.g. Pizzolato et al. 12003b . 
We determine with the diff'erential-rotation code of Ktiker & 
Stix (2001 ) as a function of depth based on a description of the 
pressure scale height. Hp, and the classical mixing length, a. The 
characteristic convective length is then defined Sis£ = a Hp, with 
a = 2. In the case of a fast rotator, the convective heat flux is no 
longer aligned with the temperature gradient and tilts towards 
the rotational axis. It then creates a horizontal heat flux from the 
equator to the poles (meridionally circulated) that, together with 
the Reynolds stress, maintains difl'erential rotation. Fig. [6] shows 
the turn-over time for HD 123351 plotted versus the fractional 
stellar radius. The depth of the convection zone is predicted from 
the MESAQ code (Paxton et Sil.WlU) to be 90% of the stellar ra- 
dius. The average turn-over time is around 40 d. It would give a 
Rossby number of ^ 1 .4 and, according to the dynamo criterium 
of Durney & Latour ( 1978 ), be insufficient to operate an eff'ec- 
tive dynamo, in contradiction with the observations. If we de- 
termine the Rossby number using the convective turnover time 
at the base of the stellar convection zone, i.e. Ro ^0.13 with 
Tc ~ 80 d, it agrees with the dynamo criterium. This may indi- 
cate that the dynamo of HD 12335 1 is indeed located at the base 
of the convection zone, as in the Sun. 

3. 7. Chromospheric activity 

We measure an average equivalent width of the Hei D3-line of 
39 mA with an estimated error of 5 mA, which is comparable to 
the measurements for most active G and K dwarfs in the sam- 
ple of Saar et al. ( 1997) . The detection of the Hei D3-line at 
587.56 nm indicates the presence of non-thermal heating pro- 
cesses on HD 123351, since the line cannot be excited at pho- 
tospheric temperatures. It is thus a comparably "pure" indicator 
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I 4«10^- 
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<D 
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<D 
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Fig. 6. Convective turn-over times, r^, in seconds as a function of 
depth in the convection zone of HD 123351. Note that the turn- 
around of the curve below a fractional radius of ^^0.35 strongly 
depends on the adopted stellar model and could be an artefact. 
Tc spans between a maximum of 80 d near the bottom of the 
convection zone and 10 d near the surface. 



Table 5. Summary of period determinations. 



Indicator 


Period 


Technique 


FAP 




fdavs) 






y/-band photometry 


61.9+0.4 


DFT 




Radial velocities 


147.8919 


Orbit 




XvCSiLlU-Cll idUidl VClUCiliCS) 


61.0+1.1 


DFT 




Spot migration 


58.32+0.01 


LCM 




Hof line-core flux 


58.7+0.2 


DFT 




Can H&K flux 


63+5 


DFT 




Can H&K flux 


^147 


DFT 




HeiD3 


60.68+0.15 


Lomb 


1 10-6 


Magnetic flux 


58.7+4.3 


Lomb 


310-6 


Magnetic flux 


144.3+8.9 


Lomb 


0.1 



DFT: discrete Fourier transform. LCM: light-curve modeling. Lomb: 
Lomb-Scargle. FAP: false alarm probability. 



of chromospheric magnetic activity because of the absence of 
photospheric contamination present in most other optical activ- 
ity diagnostics, e.g. Can H&K and Ha. However, the compara- 
ble weakness of the line makes it more likely to be affected by a 
limited S/N ratio. 

Besides Hei D3, both Can features, the H&K doublet, and 
the infrared triplet, appear as moderately- strong central emission 
lines, while Ha appears as a filled-in absorption line with an 
equivalent width of 850 mA. Absolute Ca ii H&K emission-line 
fluxes for HD 12335 1 are determined with the method of Linsky 
et al. (11979b as adapted and applied, e.g., in Strassmeier et al. 
(119901) . It is based on the measurement of the relative flux in a 
5.0-nm bandpass centered at 395.0 nm, and the relative flux in 
the two emission lines, defined by the blue and red emission-line 
footpoints, respectively. Absolute fluxes are then determined us- 
ing the Barnes-Evans relation (Barnes & Evans 1976) and the 
measured value of the star's Johnson V - R color. Following 
the recipes in Hall (1996), the relative fluxes convert into abso- 
lute chromospheric emission-line fluxes of, on average, 1.6x10^ 
ergcm"^s"^ for the sum of the Can H and K lines, 1.4x10^ 
ergcm"^s"^ for Ha and 8x10^ ergcm"^s"^ for the sum of the 
Can infrared-triplet lines. These are values typical of RS CVn 
stars and indicative of a magnetically active atmosphere. 

4. Time-series analysis 

4.1. The location and evolution of photospheric spots 

The spot longitude, which is the most stable and well-determined 
parameter, arises directly from the light-curve shape and marks 
the presumed spot concentration. The longitude difl'erence be- 
tween the spots can depict systematic surface patterns. When the 
spots are very close, then a single, elongated spotted region prob- 
ably causes the light variation that then appears quasi- sinusoidal. 
On the other hand, double-humped light curves indicate spotted 
regions on the two hemispheres of the star separated by nearly 
180°. We perform a more detailed light-curve fitting with our 
time-series spot-modelling code SpotModeL. The strengths and 
shortcomings of the code were described in detail in Ribarik et 
al. (2003 ). The program reconstructs the position, size, and tem- 
perature of up to three cool spots by minimizing the fit residuals 
with the help of the Levenberg-Marquardt algorithm as formu- 
lated in Press et al. ( 2002). Free parameters are the spot radii, 
the temperature relative to the (fixed) photospheric temperature, 
and their longitudes and latitudes. In the present application, we 
solve for the V and Ic Hght curves simultaneously. Figure[T]illus- 



trates the fits along with a graphic presentation of the resulting 
two spot models. 

We first proceed with the determination of the spot temper- 
atures and the best-fit spot latitudes. For this initial step, we as- 
sume the average photometric period (61.9 days) to be the ro- 
tational period. Average spot temperatures are obtained by a si- 
multaneous two- spot fit of the seasonal V and Ic Hght curves. 
The seasonal light curves usually span three stellar rotations with 
significant changes in the light-curve shape in V and in Ic. Spot 
latitudes are kept as a free parameter in this initial step. We note 
that these spot latitudes are strongly model dependent. Here we 
assume two circular spots. The average rms of the fits to V - Ic 
implies spot temperatures precise to roughly ± 100 K, with a span 
of 20-180 K for the various seasons. Any seasonal change be- 
low this threshold is deemed unrealistic from our data. The grand 
average spot temperature over the 12 years was 3620 K (rms of 
250 K), i.e. cooler than the unspotted photosphere by 1160 K. 
With the exception of years 2000 and 2009, all seasonal rms val- 
ues are below or around 100 K, i.e. 1.4 to 8 times smaller than 
the 250 K from the grand average, thus considered real diff'er- 
ences. In a second step, we again applied SpotModeL but now 
to the V and Ic Hght curves in the time domain and with the spot 
temperatures and spot latitudes kept fixed to the corresponding 
average seasonal values listed in Tabled In that table, £ denotes 
the spot longitude, b the spot latitude (positive in the dominating 
hemisphere), A the fractional spot coverage in per cent of the to- 
tal sphere, Tspot the spot temperature, and (Xfit the squared sum of 
the residuals in V and / in mmag. Labels 1 and 2 denote spots 1 
and 2, respectively. 

The next step included the revision of the rotational period. 
This was done by determining the slope of the longitudinal spot 
migration in the longitude versus (vs.) time diagram for each 
spot and during each observing season. Then, the seasonally av- 
eraged longitude of each spot, or spot group, was allowed to be 
shifted by an integer rotation such that the resulting long-term 
migration was brought into the best-possible agreement with the 
seasonal migration (Fig. [7^). This was achieved for a slope of 
-0.38 deg /d for spot 1 and -0.39 for spot 2, corresponding to 
a new period of 58.32 days. As a comparison, the seasonal aver- 
age slopes were -0.32 and -0.35 for spots 1 and 2, respectively. 
The slight diff'erence between the two spots is due to spot 1 ap- 
pearing at a slightly longer longitude than expected in 2000 and 
by spot 2 being at a markedly shorter longitude than expected 
in 2007. We tend to ignore the diff'erence and adopt the aver- 
age slope as the true rotation period of the star. We then applied 
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Table 4. Seasonal average spot parameters 1998-2009. 



Year 


mid HJD 


^1 




b^ 


bo 




An 


T 

-'^ spot 


^ tit 






C) 


C) 


C) 


C) 


(%) 




(K) 


fmmas) 


1998 


2450943.815 


313.8+2.0 


163.5+4.3 


56+4 


55+9 


7.45+0.08 


2.93+0.07 


3830+20 


6.0 


1999 


2451264.194 


307.9+3.0 


105.7+2.3 


52+6 


53+4 


4.07+0.06 


5.83+0.07 


3880+25 


5.6 


2000 


2451620.425 


202.0+4.1 


61.5+7.9 


55+8 


20+22 


4.42+0.07 


1.91+0.17 


3250+135 


12.4 


2001 


2452026.519 


257.6+2.3 


52.1+4.4 


61+5 


20+10 


4.31+0.03 


2.19+0.09 


3210+75 


9.3 


2002 


2452359.676 


238.5+2.7 


52.1 + 1.7 


23+8 


55+4 


1.88+0.04 


3.61+0.02 


3500+60 


7.6 


2003 


2452702.954 


223.6+1.0 


63.7+3.4 


48+3 


12+8 


3.83+0.03 


1.40+0.05 


3420+60 


5.0 


2004 


2453060.370 


249.2+2.1 


75.2+4.4 


41+5 


26+9 


3.23+0.03 


1.57+0.05 


3520+90 


6.8 


2005 


2453444.012 


167.5+2.3 


47.8+4.2 


20+5 


32+6 


6.29+0.18 


3.36+0.12 


3850+80 


7.6 


2006 


2453830.231 


183.6+1.7 


43.7+5.4 


36+4 


18+10 


4.30+0.06 


1.51+0.09 


3525+100 


6.1 


2007 


2454198.866 


184.4+2.3 


-37.1+3.0 


28+5 


44+5 


4.00+0.10 


3.05+0.07 


3600+70 


6.8 


2008 


2454567.950 


204.8+1.5 


22.5+4.1 


26+9 


66+14 


6.46+0.06 


2.86+0.03 


3970+70 


5.5 


2009 


2454885.704 


245.4+7.5 


30.7+7.7 


30+15 


42+17 


5.07+0.25 


4.72+0.19 


3870+180 


20.5 



Note. Both spots have equal temperature. The errors come from the summed residuals from the simultaneous light-curve fit in V and Ic- The 
stellar rotation period was fixed to 58.32 d. 



1999 2001 2003 2005 2007 2009 
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Fig. 7. Results from the photometric spot modeling, a The reconstructed spot longitudes expressed in stellar rotations versus time. 
Filled gray points refer to spot 1, open circles to spot 2. These values were obtained based on the initial photometric period of 61.9 d 
and after shifting full seasons by an integer rotation such that the resulting long-term migration is brought into agreement with the 
seasonal migration (the latter is indicated by the short lines). The two lines are linear regression fits and their slopes indicate the 
true rotation period of 58.32 d. b Top panel: spot longitudes versus time based on the 58.32-d period. Filled gray points refer to 
spot 1 , open circles to spot 2 as in panel a. The pluses indicate the longitudinal separation of the two spots and show a fairly constant 
value of 150° ±17°. Bottom panel: the reconstructed fractional spot area versus time. In both panels, seasonal average values are 
additionally plotted. 



SpotModeL again in the time domain but with the rotation pe- 
riod fixed to 58.32 days. This led to the refined spot parame- 
ters in Table m We emphasize that our spot solutions, although 
precise and consistent, are not a detailed accurate representation 
of the real stellar surface because of the assumption of circu- 
lar spots and the mathematically non-unique fit (see Savanov & 
Strassmeier 2008} but serve as a qualitative surface map for fur- 
ther discussions. 

Fig. [7J) shows the migration of the spot longitudes and the 
spot size together with the (longitudinal) separation of the two 
spotted regions based on a period of 58.32 days. We note that 
in all cases only two spots were needed to fit the light curves 
to within their observational errors. We interpret them as active 
longitudes that appear to migrate in both directions with respect 
to the average photometric period but with more-or-less constant 



longitudinal separation. There is a short-term rapid migration of 
individual spots (appearing as jumps in Fig.fTJ), top panel), e.g. 
of spot 1 between 1999 and 2000. This migration during the ob- 
serving season 2000 was paralleled with a systematic increase 
in its area by ^^40% as well as a decay in the area of spot 2 by 
even a factor of three from an almost equal area at the beginning 
of 2000 to almost disappearance at the end. During 2000-2004, 
the separation between the spots increased and reached a max- 
imum separation in 2004 when indeed a double-humped light 
curve was observed. The rapid migration of spot 1 around 2005 
caused the two spots to become closer and after that the sepa- 
ration increased again, reaching an equal maximum in 2008 as 
in 2004. From the time-frequency analysis in Fig. [3]3, we find 
that around 2004, i.e. when the spot longitude separation was 
the largest, the photometric period had split into two distinct val- 
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Fig. 8. Modulation of chromospheric activity tracers during the years 2007-2010. The two top rows show Ha (a), the two bottom 
rows show Can H&K (b). The left panels are overplots of all line profiles for 2007-2010 (in panel b only the K-line is shown), 
while the bright lines indicate the median profiles, respectively. The right hand panels are the dynamic spectra phased with the 
photometric rotational period (top) and the orbital period (bottom) and the zero points from Table [2l The color scale indicates the 
relative emission-line flux with respect to the median profile. The two dashed, horizontal lines show ±v sin /, and the chain of points 
in the very bottom of the graphs represents the phase coverage (one dot is one spectrum). It appears that both activity tracers display 
a clear modulation with the rotational period, with a tendency to cluster around a time of periastron at phase 1.0, and maintained 
that behavior for the entire duration of our observations. 



ues (62% and 57^3). This suggests that in reality the two spotted 
regions were at difl'erent latitudes on a difl'erentially rotating sur- 
face with AP/P ^0.076 as a lower limit, albeit a short-term, fast, 
spot-morphology change could not be excluded as an alternative 
explanation. We note that the sign of the difl'erential rotation can- 
not be reconstructed because of the inherent lack of a reliable lat- 
itude determination from one-dimensional photometry. We note 
again that the longitudinal separation of the two spots remained 
nearly constant during the 12 years of observations (150°+ 17°). 
This suggests that the two spots, or spotted regions, have a com- 
mon origin inside the star and are not just two independent flux 
tube rings. However, we see no evidence of a periodic flip flop 
of the active longitudes, as found for many other spotted stars 



with long-term photometric coverage (e.g. Korhonen & Jarvinen 
[20Q7b . Only in 1999 and in 2002 was spot 2 larger than spot 1. 

4.2. Modulation of chromospheric tracers 

Our STELLA spectra contain a number of chromospheric activ- 
ity tracers, most notably Ha and Can H&K. With a total time 
span of four consecutive years, 2006-2010, and a sampling of 
the order of one spectrum per (clear) night, the spectra are suit- 
able for searching a modulation as short as a fraction of the star's 
rotational period. Table [5] summarizes the results. 

Balmer Ha. The spectra were first re-normalized to optimize 
the continuum setting in the echelle order that contains the Ha 
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line. We computed a median spectrum from the total of 955 spec- 
tra and then iteratively shifted the individual spectra in intensity 
and wavelength by minimizing the scatter on both sides ±0.3nm 
around the Ha line. Fig. [8^ is a time series plot of the Ha-line 
profiles for the entire time series together with a dynamic plot 
of its residual spectra after removing the median line profile. We 
then determined the relative 0. 1-nm flux centered on the Ha lab- 
oratory wavelength and subjected this to a period search with a 
discrete Fourier analysis. This search revealed the strongest peak 
for a period of 58.7+0.2 d, in agreement with the 58.32-days ob- 
tained from our photometric analysis (note that DFT errors are 
internal errors). After prewhitening with this period, two more 
periods become significant, one around 90 d and another at the 
orbital period of 147 d. Because 1/147 + 1/90 ^ 1/58, we con- 
clude that the ^^90-day peak is an alias and that we see very weak 
modulation of Ha also with the orbital period. The latter is only 
marginally significant though with an amplitude of around one 
standard deviation. The two panels in Fig. [9]3 show plots of the 
line-core emission phased with the rotational period and the or- 
bital period, respectively. 

Ca II H&K. We construct relative Ca ii H&K emission-line 
fluxes based on the measurement of the flux in a 5.0-nm band- 
pass centered at 395.0 nm, F50, and a 0.1-nm flux in the 
two emission lines, Fn and Fk, respectively. The sum of the 
emission-line fluxes with respect to the 5.0-nm flux, (Fu + 
Fk)/F5o, as well as the ratios of each of the two emission lines 
are then subjected to a period search similar to Ha. Fig. (8)3 shows 
dynamic spectra of both emission-line profiles with a pattern in 
agreement with Ha, albeit a bit "cleaner". The S/N is lower by 
a factor of two relative to Ha because of the lower efficiency of 
the SES at blue wavelengths and the star being fainter in the blue 
by one magnitude. Nevertheless, the periodogram from the four- 
year data set reveals a clear peak at a period of 63+5 days, which 
is in moderate but formal agreement with the periods from Ha 
and the optical photometry. Surprisingly, the 90-day alias and 
the 147-day peaks are more evident in Can H&K than in Ha, 
although less well-defined for the combined data in 2007-2010 
than for observing season 2008 shown in Fig. [9]:. It indicates 
variations from year to year. The two panels in Fig. [9]: are phase 
plots of the relative emission-line flux (relative to the median 
profile) with the rotational and the orbital period, respectively. It 
shows a ^^l-cr drop in flux during periastron passage at around 
phase or 1 . 

He I D3. For this indicator, we chose to simply fit a Gaussian 
to its line profile and then use the Gaussian's peak residual in- 
tensity for the period search. Another approach, similar to mea- 
suring the equivalent widths directly between two defined points 
in the continuum, did not ensure in the same homogeneity of the 
results. This is mostly because the line is rather weak and narrow 
with a full- width at the continuum level of about 0.014 nm, with 
the red line wing blended by a nearby Fei line (587.630 nm). 
As we do not resolve the line width, a Gaussian fit appears 
to be an acceptable compromise. A Lomb-Scargle periodogram 
(e.g. Lomb 1976) of the Gaussian-fit residual intensities from 
all epochs combined gives a well-defined peak with a period of 
60.68+0.15 d, i.e. very close to the expected rotational period of 
the star. A weighted wavelet transform analysis (Foster 1996), 
i.e. without a re- sampling of the data along the time axis as for 
the Lomb method, shows that this period is present during all 
seasons. However, it is always 10-15(T longer than the true ro- 
tation period from spot migration. The orbital period does not 
appear as clearly as the rotational period. This is partly because 
of the shorter exposure time that had been required resolve the 
rapid change in the radial velocity during periastron introducing 
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Fig. 9. Residual variations for the observing season 2008. a V- 
band light curve, b Hof-core flux, c relative Ca 11 H&K emission- 
line flux, and d He i D3 absorption line depth. The left column 
shows the data phased with the rotational period of 58.32 d, the 
right column with the orbital period of 147.89 d, both with re- 
spect to a time of periastron passage. Units are as in Fig.[8l We 
note the sharp decline in Ca 11 H&K flux around periastron in 
2008. For better visibility we plot the phases twice. 



phase-dependent scatter. Fig. [9]l shows the residual intensity of 
the line core, again phased with the 58.32-d rotational period and 
the orbital period, respectively. 
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4.3. Magnetic analysis 

A single STELLA/SES spectrum provides a mean S/N of about 
125, too low for a direct magnetic analysis. However, this can 
be compensated for by our having 273 well- sampled spectra 
in 2008 alone. This time sampling of each individual spectral 
line profile can be used in a principal component analysis (PCA; 
Bishop 1995 ) approach to boost the S/N of a particular line by 
a factor of more than five. It also allows us to search for phase- 
dependent variations in the Zeeman- sensitive spectral line pro- 
files. To diff'erentiate the small excess line-broadening caused by 
a surface magnetic field from other line-broadening mechanisms 
such as thermal eff'ects, micro- and macroturbulence or atomic- 
line parameter variation, e.g. of the oscillator strengths, we need 
to use a sophisticated line-profile inversion technique. Examples 
of this method that is based on radiative transfer in magnetic at- 
mospheres has for a long time been applied to various active cool 
stars (e.g., Saar et 1986, Valenti et al. .l995, Ruedi et al. lT997l 
Anderson et al. 120101) . 

4.3.1 . Signal-to-noise ratio enhancement 

A straightforward statistical method that can be readily applied 
to individual line profiles is the technique of PCA. The appli- 
cation of PCA in the present work is twofold. Firstly, it pro- 
vides a simple means of detecting the dominant variations in the 
data set by building the sample covariance matrix of the time- 
sampled measurements for each individual spectral line profile 
of interest. This is done in the following separately for the two 
spectral lines Fei 617.3 nm and Feli 624.0 nm. Calculating its 
eigenvalues and eigenprofiles allows us to capture the dominant 
(energy) variances in the leading eigenprofiles (principal compo- 
nents). The second motivation for using PCA is its noise-filtering 
capability. By retaining only the leading components in the data, 
one can obtain a significant reduction in noise without compro- 
mising the basic physical content in the line profile. Details of 
the application procedure to Stokes profiles were described in 
Carroll et al. f^Of) and Martinez et al. ( 20080 . We note that 
the PCA analysis in this work diff'ers from these studies in that 
we build the covariance matrix over the time domain, i.e. for a 
particular wavelength (target line) many line profiles measured 
over time are used rather than many spectral lines with diff'er- 
ent wavelength at a certain time. Applying PCA to our STELLA 
spectra and reconstructing the individual line profiles with just a 
small number of principal components (three eigenprofiles), we 
reached a formal S/N of 650, i.e. an enhancement of approxi- 
mately a factor of five with respect to the original data. 

4.3.2. Magnetic flux density inversion 

Our aim is to obtain the magnetic-flux density over the respective 
visible surface of the star. The magnetic flux density, or loosely 
called magnetic flux, is described by the product B.f, the mag- 
netic field strength B times the surface filling factor /. In the 
present work no attempt is made to distinguish these two quan- 
tities, our main goal being to provide an estimate of the mag- 
netic flux and its time variation. All inversions are performed 
under the common assumption that the magnetic field is homo- 
geneously distributed over the surface and radially aligned. 

We follow the analysis of Ruedi et al. (119971 ) and choose 
the Zeeman- sensitive line Fei 617.3 (eff'ective Lande factor 
^eff=2.50) and the Zeeman-insensitive line Fe 1 624.0 (geff=1.00) 
for a simultaneous analysis of the two lines. This iron line pair 
has the advantage of being of similar strength and having the 
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Fig. 10. Inversion fits to the averaged line profiles for the 
observing season 2008. Panel a shows the Zeeman- sensitive 
Fei 6173 line (eff'ective Lande factor geff=2.50), panel b shows 
the Zeeman-insensitive Fei 6240 line (geff=1.00). The dots are 
the PCA-enhanced observations, the lines are the inversion fits 
with a total of six free parameters. The respective lower pan- 
els show the residuals enhanced by a factor of ten. The average 
magnetic flux density is 542+72 G. 



same excitation potential 0^^=2.22 eV) so that their response 
to temperature variations can be expected to be similar. Owing 
to the temperature inhomogeneities across the stellar surface of 
HD 12335 1 , a simultaneous inversion of this line pair is expected 
to reduce the cross-talk between temperature-related broadening 
and magnetic-field eff'ects. All spectral line blends found in the 
VALD database (Piskunov et al. 1995 ) with a line depth of more 
than 0.02 within a 0.03 nm range about the line center of both 
iron lines are taken into account. 

The flux inversion is done with a variant of our Zeeman- 
Doppler imaging code iMap (Carroll et al. 2007). The forward 
module of this code consists of an LTE polarized radiative- 
transfer computation with the DELO integration method and is 
described in more detail in Carroll et al. ( 2008 ) to which we refer 
for more details and references. In addition to the magnetic flux 
density B.f, we consider the efl'ective temperature of the atmo- 
spheric model as a free parameter. The atmospheres are Kurucz 
(1993 ) ATLAS-9 models. Spectral lines are calculated in the L- 
S coupling regime and atomic line parameters are taken from 
VALD, except for the logarithmic oscillator strength, which is 
determined by an initial inversion fit of the average spectral line 
profile. Assuming the oscillator strength as a free parameter fol- 
lows the argument of Basri & Marcy (1988) and Ruedi et al. 
(1997 ) that solar values are inappropriate for the analysis of dif- 
ferent stellar types. The optimization procedure relies again on 
the Levenberg-Marquardt algorithm. 
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The underlying stellar model used in the inversion consists of 
a segmented surface with a 5° x 5° resolution. Each segment ob- 
tains the magnetic flux density and the efl'ective temperature of 
the model atmosphere. During the disk-integration process, the 
atmospheric model with its depth stratification of the tempera- 
ture and pressure is adjusted according to the line-of-sight angle 
to account for limb-darkening in the flux profiles. In total, the 
inversion solves for the following six free parameters for both 
spectral lines simultaneously; the magnetic flux density B.f, the 
eff'ective temperature, the microturbulence, the radial-tangential 
macroturbulence with equal radial and tangential components, 
an instrumental broadening modeled with a Gaussian velocity 
distribution, and the oscillator strength. The rotational velocity 
and the iron abundance are taken from the pre-analysis in Sect. [5] 
as they are not expected to be modulated. 

4.3.3. Results from the inversion 

We concentrate on the data set obtained in the year 2008, as was 
the most finely sampled. Selecting only the highest S/N spectra 
still gives a database of 273 spectral line profiles for each of the 
above iron lines. In the first step, we only invert the average line 
profiles (zeroth eigenprofile) that determine the global surface- 
averaged values of B.f and the other free parameters. This gives 
a first hint on the average strength of this field. In addition, we 
can determine the other free parameters and then keep them fixed 
for the following analysis of the individual line profiles. 

From Fig.[TOl we see that the fits to the averaged line profiles 
are extremely precise and consistent with the low noise level and 
a reduced error of 1 . 12 for the Fe 1 617.3 line, and 1.58 for the 
Fei 624.0 line. The obtained (logarithmic) oscillator strengths 
are -2.78 for the Fe 1-617.3 line, and -3.18 for the Fe 1-624.0 line. 
Micro- and macroturbulence were found to be 1.15 kms"^ and 
3.9 km s"\ respectively, and the average eff'ective temperature to 
4830 K (with a peak- to- valley variation of 60 K). The inferred 
magnetic flux density was 542 G with a formal error of 72 G cal- 
culated from the covariance matrix (see Carroll et al. 12007b . We 
emphasize that the continuum depression caused by a slightly 
reduced efl'ective temperature related in turn to a non-zero spot- 
filling factor, is implicitly solved for in our inversion. If we plot 
the eff'ective temperature of each line pair versus rotational (or 
orbital) phase, it shows plots similar to those in Figs. [TTh and 
b. This is expected if the magnetic flux density stems from lo- 
calized regions rather than from a global network. An indication 
of the reality of the magnetic field broadening is obtained from 
a number of inversions without a magnetic field. In all of these 
zero-field cases, the quality of the fits ix^ of 1 .74 and 1 .92 for the 
two lines, respectively) were unable to reach that of the magnetic 
case. This is also obvious in the projected landscape, which 
shows that the zero-field case is far beyond the 5 -a region of the 

minimum with a magnetic field. Therefore, we consider our 
magnetic-flux result as significant. 

The next step in our analysis is to invert the 273 individual 
line profiles of our iron-line pair to search for a systematic peri- 
odic behavior. Fig. [TT] shows the results and indicates a clear ro- 
tational dependence on the magnetic flux, as well as a (weaker) 
dependence on the orbital phase. The individual formal errors in 
the retrieved magnetic flux values are between 55 and 80 G, de- 
pending on the initial S/N. We note again that by performing a 
simultaneous inversion of the two iron lines with both having a 
similar temperature response but a difl'erent magnetic response, 
we reduce the crosstalk of temperature efl'ects and its possible 
misinterpretation in the inversion process. A Lomb-Scargle pe- 
riodogram (Fig.fTTb) indeed indicates a pronounced peak above 



a 0.01% false-alarm probability at 58.7 d with an uncertainty 
of 4.3 d determined from a Monte Carlo resampling. A second 
significant peak, but not quite with as high a probability as the 
first peak, shows a periodicity of 144.3 d with an uncertainty of 
8.9 d. The distribution of the summed squared errors for each 
line-profile fit versus phase does not indicate any systematic dis- 
tribution, suggesting that the inversion was not fooled by a sys- 
tematic, residual surface-temperature variation. As in Fig. [9] for 
the chromospheric indicators, we plot the magnetic flux densi- 
ties over the rotational and orbital phase in Fig.fTThb. 

5. Discussion 

5.1. Circularization and synchronization 

Tidal interaction between stars can couple their spin to their or- 
bital angular momentum such that dissipation of tidal energy 
causes a binary to circulari ze its o rbit (e.g. Zahn 1989 and earlier 
papers; Tassoul & Tassoul [19961). Observational results indicate 
that most active binaries have circular orbits with synchronized 
rotation and a spin a xis perpendicular to the orbital plane (e.g. 
Fekel & Fitter 119891 ). Other systems with weaker tidal interac- 
tions may not be currently old enough to have reached complete 
synchronization or circularization. However, these conditions 
may occur once tidal interactions have operated long enough or 
after a system evolved to a state in which tidal interactions are 
enhanced. For late-type stars with convective envelopes, Zahn 
(1989) investigated the efl'ects of the equilibrium tide on syn- 
chronization and circularization, while Tassoul (119881) explored 
the theory of geometrical distortions that cause large-scale hy- 
drodynamic currents. Although these two theories disagree sig- 
nificantly on absolute timescales, both predict that synchroniza- 
tion should occur before circularization. 

HD 12335 1 is neither synchronized nor circularized. Its large 
orbital eccentricity suggests a pseudo-synchronous rotation pe- 
riod of 12 d according to the tidal-friction theory of Hut ( 1 9?T]) , 
which is just 8% of the orbital period. With a 58. 3-d rotation pe- 
riod, the star appears to be a strongly asynchronous and slow 
rotator. Since all convective stars likely undergo some sort of 
diff'erential surface rotation, one would expect a certain range 
of rotation periods if derived from photometric spot activity be- 
cause spots could appear at various stellar latitudes. For com- 
parison, the (strong) solar diff'erential rotation allows a spread 
of rotation periods of ^10%. Comparable results are found for 
other active stars based on Doppler imaging (e.g. Barnes et al. 
12005 1 Strassmeier 2009). However, the diff'erence between the 
observed photometric period and the pseudo-synchronous rota- 
tion period for HD 123351 is more than 500% and could not be 
explained by diff'erential surface rotation. 

Therefore, we ask what caused a cool star to fail to reach 
equilibrium in terms of its spin and orbital angular momentum 
after 6-7 Gyr? Although not in disagreement with the predic- 
tion that synchronization shall occur before circularization, the 
Tassoul and Zahn mechanisms predict vastly diff'erent timescales 
for HD 123351. While the original Zahn (fT989b equations pre- 
dict synchronization and circularization well beyond the age of 
the universe, the Tassoul (1988) equations predict synchroniza- 
tion after ^^200,000 years and circularization after ^^3 Gyr. Both 
can obviously be excluded because of the measurements derived 
from the observations. Additional mechanisms were proposed 
to account for the strong tidal coupling required by the obser- 
vations of circularization of binary systems in open clusters of 
various age (Ogilvie & Lin 2007b , but we refer to the refutation 
of Tassoul 's hydrodynamical mechanism by the proper account- 
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Fig. 11. Magnetic flux variation for the observing season 2008. a Phase plot with the stellar rotation period of 58.32 d, b phase 
plot with the orbital period of 147.89 d, and c Lomb-Scargle periodogram from the 273 individual spectral line inversions. The two 
strongest peaks in the periodogram agree with the rotational and the orbital periods, respectively. The strongest peak refers to a 
period of 58.7 d, the second strongest peak to 144.3 d. Note again that we plot the phases twice for clearer visibility. 



ing of an Ekman layer by Rieutord & Zahn (119971) . In the case 
of HD 1 2335 1 , the orbital and rotational frequencies today imply 
that inertial oscillations (waves) might have neen excited in the 
star leading to an average tidal dissipation approximately 10- 
100 times that predicted by Zahn's theory. Such inertial waves 
are not described in the traditional approximations and corre- 
spondingly reduce the estimated synchronization and circular- 
ization times. Because non-linearities are expected for eccentric 
binaries and for compone nts old er than the Sun, no simple pre- 
diction of Ogilvie & Lin ( 2QQ7b is possible for HD 123351, but 
in general brings the theory closer to the observation. 

5.2. Inter-binary activity? 

The assumption of an inter-binary magnetic field that could 
transfer orbital to spin momentum and vice versa may indeed 
be inaccurate for intervals longer than purely gravitationally and 
hydrodynamically predicted timescales. It is well known that 
inter-binary magnetic fields focus the accretion stream in cata- 
clysmic variables (CV) containing a magnetic white dwarf (e.g. 
Warner & Woudt 120021) , as well as some Algol systems (e.g. 
Richards i2007jj). Reconnection processes within these fields may 
also act as the source of inter-binary flares, suggested e.g. for 
the RS CVn binary system HR 5110 (Grafl'agnino et al. 17995 ). 
The field lines in general provide a channel for the exchange of 
charged particles between the two stellar components. Its princi- 
pal existence on a geometrically shorter scale is nicely demon- 
strated by the Jupiter-Io system (e.g. Su 2009). However, the 
closest distance between the two components in the HD 123351 
binary system (^^5 R^) is still larger than for the aforementioned 
CV and Algol binaries, thus the expected field density may be 
insufficient for inter-binary flares. No such superflare was de- 
tected in our data. Coronal mass ejections may however still ex- 
ist, which would have had an effect on the evolution of the spin 
and orbital angular momentum of the star but remains specula- 
tive. 

5.3. Orbital-induced magnetic activity? 

In a very eccentric close binary similar to HD 123351, we might 
have the rare possibility of exploring a hypothetical efficiency 
increase and decrease of a stellar dynamo when the two com- 
ponents are closest and farthest, respectively, assuming e.g. that 
some sort of extra convective mixing coupl es the dynam o ac- 
tion to tidal forces. Glebocki & Stawikowski fl977||T288| found 
some relation between chromospheric Can emissions and the 



orbital parameters of late-type giants and suggested that tidal 
forces are somehow responsible for the enhanced chromospheric 
activity. In contrast, Basri (1987) ruled out tidal coupling as an 
explanation of the increased activity in binaries. Schrijver & 
Zwaan (119911) argued that in a synchronized close binary the rel- 
evant axis for rotational effects is through the center of gravity 
of the binary system and thereby could affect the stellar dynamo 
in such a way that it produces enhanced magnetic activity rela- 
tive to single stars. Zaqarashvili et al. (2002) assumed that the 
deviation from a spherical shape of a binary component due to 
tidal interaction excites a fundamental pulsation mode that then 
amplifies the torsional waves in the interior magnetic field and 
eventually causes greater surface activity. It has been found that 
(extra-solar) planets can induce chromospheric and coronal ac- 
tivity in the form of a hot spot (Shkolnik et al. 12003b , although 
not all stars with close-in planets are prone to these star-planet 
interactions (e.g. Fares et al. 2010). Lanza (2 008 , 2009) argued 
that a flux rope between the host star's coronal magnetic field 
and the planet's magnetosphere experiences a decrease in rel- 
ative helicity and reacts with an energy release from the flux 
rope's footpoint(s). His model predicts a greater prominence ac- 
tivity for stars with planets than without, caused by the accumu- 
lation of matter evaporated from the planet. A very similar pic- 
ture may be adopted for binaries such as HD 123351 that have 
a low-mass dwarf star as a companion instead of a close-in hot 
Jupiter. 

Our Can H&K and Hei data indeed indicate a correspon- 
dence of magnetic activity with orbital periastron, at least dur- 
ing observing season 2008. No clear relation was seen for the 
observing seasons 2007 and 2009. However, this may be caused 
by some bias because the season 2008 was the one with the 
highest data quality and finest data sampling, both of which are 
needed to detect such a correspondence. Furthermore, the mod- 
ulation amplitudes were always small and on the order of the 
annual rms. But the detection is quite surprising because Ha 
does not show a similar correspondence at all, which is the op- 
posite of what we would expect if there is prominence activity 
(e.g. Mackay et al. 2010, Frasca et al. 2008 ). Our magnetic-flux 
analysis also revealed evidence of a weak orbital-phase coher- 
ence in 2008 and thus strengthens the Can and Hei detections. 
The B.f orbital variation is in phase with the variation in both 
Can and Hei in the sense that lower chromospheric emission 
is seen when the (photospheric) magnetic flux density is lower, 
with their respective minima coinciding with a time when the 
two binary components are closest in distance, i.e. at periastron. 
The orbitally phased V-band light curve in Fig. [9^ is completely 
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out of phase, as expected because it is purely photospheric in 
origin. The comparison of the chromo spheric emissions with the 
magnetic flux density is more puzzling when plotted with the 
rotational period. Then, B.f varies exactly in anti-phase with 
all three chromospheric lines, in particular also with Ha (com- 
pare the left columns in Fig.[9]3cd with Fig.fTTh). which is unex- 
pected. The rotationally phased V-light curve appears nicely in 
phase with B.f and in anti phase with the chromospheric lines, 
as expected if chromospheric plages are superposed on photo- 
spheric spots as in the Sun. In summary, we take the anti-phase 
relation of the magnetic flux with the chromospheric emissions 
as evidence that there are two magnetic fields present at the same 
time, a localized surface magnetic field associated with spots 
and a global field that is oriented towards the secondary compo- 
nent such as an inter-binary field. We have no already developed 
model to combine all of this but conclude that the ^-component 
that we measure in B.f possibly stems from a mixture of bright 
plage-like regions and dark spots, both with diff'erent respective 
filling factors (/), but the former also caused by the interaction 
with the secondary star. 

The overall amplitude of the radial-velocity residuals shown 
in Fig.[2]3 suggests that the star was more asymmetrically spot- 
ted in 2007 than in 2008, albeit the V-light curves showed ap- 
proximately the same amplitudes and shapes. In 2007, the Ha^- 
core flux indeed showed the most pronounced variations {^Acr) 
of all seasons when phased with the rotation period but again 
basically just a scatter plot when phased with the orbital period 
(not shown). This indicates that Ha responds to the same sur- 
face inhomogeneities as seen in the other lines or the photome- 
try. It encourages us to believe that the HcK-profile variations of 
HD 12335 1 originate predominantly from near the stellar surface 
and have no detectable circumstellar (absorption) component, as 
would be expected for coronal prominence activity. We note that 
we did not apply a correction for the photospheric contribution 
to Hof but did all time-series analysis from residual line profiles 
after subtraction of a median profile. 

6. Summary and conclusions 

We have obtained a four-year-long, phase-resolved, spectral 
database of the active star HD 123351 using our new robotic 
telescope facility STELLA in Tenerife. These data are comple- 
mented with 12 years of photometric monitoring with our auto- 
matic photometric telescope Amadeus in Arizona. Our summary 
and conclusions are as follows: 

- We have discovered HD 123351 to be a single-lined spec- 
troscopic binary with a period of 147.8919 d and an orbital 
eccentricity of 0.809. 

- The STELLA Echelle Spectrograph enabled the most precise 
(and presumably also most accurate) orbit determination of 
any active binary so far in the literature with an rms of the 
residuals of just 47 m s"^ . 

- A period of ^^60 d was obtained from 12 years of Amadeus 
y/-band photometry and has been converted into a true rota- 
tion period of 58.32 d from a linear fit to the two dominant 
spot-migration curves. The period diff'erence from the two 
spots in 2004 suggests a lower limit to the surface diff'eren- 
tial rotation of AP/P of 0.076, which is roughly half the solar 
value. 

- The global rotation period of HD 12335 1 was also recovered 
from the orbital radial- velocity residuals. 

- A spectrum synthesis yielded average stellar atmospheric 
parameters of Teff =4780+70 K, logg=3.25±0.30, and 



[Fe/H]=0.00 ± 0.08. The Hipparcos parallax in combination 
with evolutionary tracks suggests a mass of 1.2+0.1 M© and 
an age of 6-7 Gyr, consistent with the measured ^^C/^^C ratio 
of ^40. 

- The lithium abundance of log n= 1 .70 appears higher than ex- 
pected for that age but the line also appears broader than ex- 
pected for the vsin/ of 1.8+0.7 kms"^ We suspect that its 
asymmetry is partly due to a significant contribution from 
the ^Li isotope. 

- From an application of the Kiiker-Stix diff'erential-rotation 
code to a MESA stellar model with matching parameters of 
HD 12335 1 , we obtain a convective turn-over time that spans 
between a maximum of 80 d near the bottom of the convec- 
tion zone to 10 d near the surface. This implies an average 
Rossby number greater than unity. 

- Our light-curve analysis has reconstructed two persistent 
spotted regions with a nearly constant longitudinal separa- 
tion of 150° +17° over the course of 12 years. The spot tem- 
peratures were on average 1160 K lower than the eff'ective 
temperature covering between 5-10% of the entire surface. 
The photometry showed no signs of a flip flop nor of flares. 

- Chromospheric emission-line fluxes of ^^10^ ergcm"^s"^ 
were measured for Can H&K, Hck, and the Can infrared 
triplet. The detection of the Hei D3-line at 587.56 nm also 
indicated the presence of non-thermal heating processes. 

- Time-series analysis of above chromospheric activity indi- 
cators yielded a clear modulation with the stellar rotation 
period but also provided evidence of a modulation with the 
orbital period, albeit less well-defined. 

- We measured the Zeeman broadening from a pair of iron 
lines and obtained an average surface magnetic flux density 
of 542+72 G. A time-series analysis of the data from 2008 
yielded again a clear modulation with the rotation period and 
a weak modulation with the orbital period in the sense that a 
lower magnetic flux was seen at times near periastron. 

- It remains unclear why the magnetic flux varies exactly in 
anti-phase with the chromospheric emissions if plotted ver- 
sus the rotational period, and is lowest near periastron when 
plotted versus the orbital period. The measured magnetic 
flux possibly stems not only from cool photospheric spots 
but also from bright plage-like regions. 

- We postulate an inter-binary magnetic field and tentatively 
conclude that it is responsible for the flux dilution at perias- 
tron. It is also possibly responsible for the unexpected slow 
and asynchronous rotation of the primary star. 
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